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Abstract 

We construct black hole solutions to the Yang-Mills equations in an A(iS'4-Schwarzschild 
background which exhibit superconductivity. What makes these backgrounds p-wave su- 
perconductors is that the order parameter is a vector, and the conductivities are strongly 
anisotropic in a manner that is suggestive of a gap with nodes. The low-lying excitations of 
the normal state have a relaxation time which grows rapidly as the temperature decreases, 
consistent with the absence of impurity scattering. A numerical exploration of quasinor- 
mal modes close to the transition temperature suggests that p-wave backgrounds are stable 
against perturbations analogous to turning on a.p+ip gap, whereas p+zp-wave configurations 
are unstable against turning into pure p-wave backgrounds. 
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1 Introduction 

In [H [2], [3l H] , three variants of a mechanism were proposed through which black holes can 
spontaneously break an abelian gauge symmetry at finite temperature, leading to a form of 
superconductivity. The symmetry breaking occurs through the formation of a superconduct- 
ing condensate that floats above the horizon. To make the mechanism work, one needs a 
charged matter field whose quanta form the superconducting layer, some interaction which 
keeps it from falling into the black hole, and some feature which prevents it from escaping 
to infinity. In the variants considered so far, the charged field is bosonic. Charged fermions 
might also work, but then one has to establish some pairing mechanism before condensation 
can occur. 

The purpose of this paper is to provide another example of superconducting black holes, 
closely related to the one in [1], and to argue that it exhibits a p-wave gap. All our results 
will be based on classical solutions to field equations of Einstein- Yang-Mills theory with a 
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negative cosmological constant: 

' ' ■ ' + , (1) 

where F":^ is the field strength of an SU (2) gauge fie 



subgroup as the gauge group of electromagnetism|l 



d. As in the plan is to regard a U{1) 
and to persuade the off-diagonal gauge 
bosons, charged under this U{1), to condense outside the horizon. Our conventions on metric 
signature and field normalizations are the same as in [Ij. The action ([T]) is almost completely 
dictated at the two-derivative level by local diffeomorphism symmetry and gauge symmetry, 
and it can be embedded in M-theoryo This endows its solutions with an interest independent 
of potential applications to superconductivity. Indeed, the solutions to be considered are 
loosely related to those of [21E], from which a significant literature has sprung, reviewed for 
example in [3, E] • 

Section [2] is devoted to a conceptual overview of black holes which superconduct through 
the mechanism described above. In section[3]we describe the background solutions of interest. 
In section m we study the electromagnetic response, along the lines of We find a frequency- 
dependent conductivity which depends strongly on the polarization of the applied electric 
field. The low-frequency behavior is suggestive of quasi-particle excitations whose dissipative 
mechanisms are entirely due to finite-temperature effects. 

In section we provide numerical evidence that the p-wave backgrounds are stable against 
small perturbations that turn on a. p + ip gap. In section \572\ we provide numerical evidence 
that the p + zp-wave backgrounds of [2] are unstable against small perturbations that turn 
them into the p-wave backgrounds described in section [31 Our numerical explorations are far 
from covering the full range of parameters, but the simplest scenario consistent with them is 
that p + zp-wave backgrounds are always unstable, and that p-wave backgrounds represent 
the thermodynamically preferred phase for T less than a critical temperature Tc. 



^More precisely, the boundary theory has a global SU{2) symmetry, and adding electromagnetism means 
weakly gauging this U{1) in the boundary theory. By contrast, the gauging of the full SU{2) symmetry in 
the gravity theory encodes aspects of the SU (2) current algebra dynamics in the boundary theory. 

^Embedding ([T]) in string theory or M-theory involves a specific choice of the gauge coupling g, and it has 
not been shown that any of the three variants [2 [U 3] of superconducting black holes occurs in M-theory. 
Also, M-theory constructions usually involve scalars that couple non-renormalizably to the gauge field and 
may be an important ingredient in constructing symmetry breaking solutions. 
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2 Conceptual overview 



In the anti-de Sitter examples of [21 [3l H] , the superconducting layer floats above the horizon 
because the horizon is also charged. Electrostatic repulsion overcomes the gravitational 
attraction that ordinarily would suck the superconducting layer into the horizon. If the 
spacetime were asymptotically fiat, then (barring some special interactions such as considered 
in [1]) one expects that electrostatic repulsion would cause the superconducting layer to be 
blown off to infinity. But asymptotically anti-de Sitter geometries prevent this. Massive 
particles, no matter how strongly repelled from a horizon, cannot reach the boundary of 
anti-de Sitter space. So they instead condense near the horizon, where "near" means that 
the field profile is normalizable, carrying finite charge if the horizon is finite, or finite charge 
density if the horizon has infinite extent. An analogy with the classic two-fiuid model of 
superconductors is possible: the charged horizon describes the normal component, and the 
condensate above it is the superconducting component. See figure [H In this analogy, it is 
important to recall that in the gauge-string duality [9], [TOl [11] , the extra dimension r is not an 
additional fiat dimension transverse to the sample; instead, it is a way of organizing energy 
scales in the dual field theory, which is strictly 2+1-dimensional and non-gravitational. Thus, 
although the condensate is "above" the horizon in the gravity picture, it interpenetrates the 
normal state in the field theory picture. 

In where the bulk geometry is based on Einstein- Yang-Mills theory, it was shown 
that there is a second order transition, with mean field theory exponents, between a non- 
superconducting state at high temperatures, where all the charge is in the normal component, 
and a superconducting state at low temperatures. A similar transition probably occurs for 
black hole geometries based on gravity coupled to the Abelian Higgs model [3l H], but this 
has not been demonstrated conclusively. What is clear, both for Einstein- Yang-Mills and for 
the Abelian Higgs model, is that the moduli space of black hole solutions includes Reissner- 
Nordstrom anti-de Sitter black holes (hereafter RNAdS), which describe the normal state jf| 
joining continuously onto a branch of symmetry-breaking solutions. The simplest argument 
supporting this picture is based on studying linear perturbations of the charged field around 
an RNAdS solution. They obey an equation of the form 



■^It has been suggested [121 [13] that RNAdS black holes are dual to a close analog of the pseudogap 
state of high Tc materials. In the context of our constructions, this does not seem quite right, because the 
fraction of charge in the condensate goes to zero near scaling as Tc — T, whereas the transition from 
superconductivity to the pseudogap state appears to take place while this fraction is finite and non-zero. 
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Figure 1: A superconducting condensate floats above a black hole horizon because of a 
balance of gravitational and electrostatic forces. The condensate carries a flnite fraction of 
the total charge density, so there is more electric flux above the condensate than there is 
right at the horizon. A massive charged particle, labeled if)^, may be driven upward by the 
electrostatic force, but because of the warped geometry of AdS^, its trajectory cannot reach 
the boundary. So ipj^ must participate in the condensate if it doesn't fall into the horizon. 
The frequency- dependent conductivity can be found by calculating an on-shell amplitude for 
a photon propagating straight down into the geometry. 



4 



where □ is an appropriate covariant wave operator and 



m 



eflf 



Here q is the charge of a quantum of the charged bosonic field ip: in the Abehan Higgs case, 
it is a complex scalar, while in the Einstein- Yang-Mills case, it is a complex combination of 
non-abelian gauge bosons. $ is the electrostatic potential, which vanishes at the horizon but 
grows quickly outside it if the electric field is strong. The metric component (7** is negative 
in the conventions we use, and it diverges to —00 at the horizon, so ([3]) implies that ip is 
tachyonic near the horizon if q is big enough and m is small enough, provided also that 
the horizon carries sufficient chargejj It is a matter of calculation to determine when ([2]) 
admits a static solution. When it does, one may reasonably assume that it signifies the 
joining of a branch of symmetry breaking solutions onto the RNAdS solutions. And one 
may calculate a critical temperature where the joining occurs. It does not necessarily 
follow that Tc is the temperature of a second order phase transition: it could be that the 
solutions which only slightly break the symmetry are thermodynamically disfavored, and 
that a first order transition to solutions with finite symmetry breaking occurs at a different 
temperature. Explicit calculations of the free energy, as in are necessary in order to 
determine the phase diagram. But even the calculations of [1] are not enough: one must 
also ask whether a solution is stable under small perturbations. At least for a certain range 
of parameters, we will show in section 15.21 that the solutions of [Ij are unstable against a 
perturbation that seems likely to turn them into p-wave solutions of the form described in 
section [31 We have not yet found an unstable perturbation of the p-wave solutions, or of the 
s-wave solutions described in [2], [3]. 

In [3] , the response of superconducting black holes to electromagnetic probes was studied. 
The black holes in question were constructed along the lines of the proposal of [2], by coupling 
the Abelian Higgs model to gravity, but a simplification was achieved by assuming that q is 
large. Provided T is not too small relative to the charge density of the horizon, this implies 
that the back-reaction of both the gauge field and the charged scalar on the metric can be 
neglected. A similar limit was considered in [1] for the Einstein- Yang-Mills case, where by 
taking the gauge coupling large and avoiding the regime of very small temperatures, one 
may similarly neglect back-reaction of the matter fields on the metric. A related limit of 



"'Actually, the most commonly considered cases have < in the case of scalars, or m = for non- 
abelian gauge bosons. The argument about massive particles' trajectories never reaching the boundary of 
anti-de Sitter space then no longer holds up, but it is replaced by standard notions of boundary conditions 
in anti-de Sitter space which again lead to normalizable condensates. 
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the proposal of [2] was studied in [T3] in the presence of magnetic field, leading to partial 
localization of the condensate in one of the spatial directions on the boundary. Also in the 
no-back-reaction limit, the effects of critical magnetic fields and vortices have been studied 
in [I51[T6]. 

The quantity of primary interest in understanding the electromagnetic response is the 
conductivity, 

cr,,M = -Gg(a;,0), (4) 

where 

GlSuj.k) = -tj d3a;e"-^^e(t)([J^(t,f), J„(0,0)]) (5) 

is the retarded Green's function of the electromagnetic current JmE The angle brackets in 
dS]) denote expectation values at finite temperature, namely 

{A) = ^tT e-^^'A Z = tT e-^^ (6) 

for any operator A. The hermitian part of aij is dissipative, while the anti-hermitian part 
is reactive^ According to a spectral decomposition, the hermitian part of Oij should be 
positive semi-definite. To see this, first note that the spacetime dependence of the hermitian 
operators Ji{t,x) is found through 

Ji{t, x) = e'^*-'^-"Vi(0, o)e-'^*+'^-^ . (7) 

Introducing a complete set of states between the two operators in ([5]) and integrating over t 
and X one obtains 

0) = - ^(2.)'.<-'(J^„)Jr"^" ^ „ . (8) 

n,m 

where 

Jr = {n\MO,0)\m) Pnm = Pn-Pm E^m = En - , (9) 

P„ and En being the eigenvalues of P and H in the state |n). Plugging ([8]) into (j4]) one 



^Our convention for indices is that i and j run over the spatial directions x and y; m and n run over 
the boundary directions t, x, and y; and fi and ly run over the bulk directions t, x, y, and r. We will later 
introduce adjoint indices a and b for the non-abelian gauge group SU{2). 

^In the theory of AC circuits it is standard to consider the complex power S = J (fxE*ji — 
J (Px E*<TijEj, whose real and imaginary parts are the real and reactive powers, respectively. The real power 
P can therefore be expressed in terms of the hermitian part of aij through P = J cPx E* ^{aij + <7*^)Ej, 
while the reactive power Q = J (Px E* -^{aij — (T*^)Ej corresponds to the anti-hermitian part. 
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straightforwardly obtains 
Z ■' ^ — ' ^ Z UJ 

n,m 

(10) 

Formally, Anm > 0, so multiplying ffTOj) by an arbitrary column vector vj to the right and 
by its adjoint v* to the left yields 



:iii 



proving that, indeed, the hermitian part of aij is positive semi-definite. 

The conductivity crjj(u;) characterizes the response to light of frequency u which is inci- 
dent on the superconductor in a direction normal to the R^'^ that the sample occupies. So 
it is perhaps intuitive that to calculate (Jij{uj) for the black hole, one should send photons 
down from the boundary of AdS^ and inquire how they are absorbed or reflected by the 
condensate and the horizon. More precisely, one uses the gauge-string duality to extract 
two-point functions from tree-level propagation of photons. The prescription for comput- 
ing such Green's functions was first articulated in [TOl [TT] . An adaptation of it to thermal 
backgrounds was correctly guessed in [T7] and then derived from the original prescription 
of [TUl [TT] in [TS] using Schwinger-Keldysh contours|l| In the case of two-point functions, 
the gauge-string prescription is closely related to D-brane black hole absorption amplitudes 
computed in a long series of papers beginning with [20]. If one expresses an asymptotically 
^^5*4 background as 

ds"^ = —i—dt^ + dx^ + dy"^) H — -dr"^ + (corrections) , (12) 

where the terms shown explicitly are the leading large r behavior, then a complexified photon 
perturbation polarized in the x direction can be expanded for large r as 



Ax — e 



(13) 



^Modulo some issues related to behavior near oj ^ 0, this prescription can also be justified [19] by the 
fact that if one analytically continues G^{uj) to the upper half-plane, then at the Matsubara frequencies 
LUn ~ 2'KnT with n > it agrees with the corresponding Fourier mode of the Euclidean correlator computed 
from the prescription proposed in |10[ 111] . 
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and the retarded Green's function is given simply by 

G«K0) = 4^, (14) 

where k = a/SttGJv is the gravitational coupling, and it is assumed that the photon wave- 
function is purely infalling at the horizon. More sophisticated examples have been discussed, 
for example, in [2T] . 

In the superconducting phase of the black holes constructed using the Abelian Higgs 
model, = o'yy and a^^y = because the order parameter is a scalar, breaking gauge 
invariance but not rotational invariance. There is a delta-function spike in Re ^^.^.(u;) at 
a; = 0, and an associated pole in Imaxxi^^)- For non-zero u and T not too close to Tc, 
Re(T^a;(u;) is very small up to a finite frequency, which can be denoted Ug = 2 A in order to 
evoke a comparison with BCS theory: A is then to be compared with the quantity denoted 
by the same letter in BCS, whose physical interpretation is the minimal energy of a single 
normal-component quasi-particle excitation. Above Ug, R.eaxx{^^) rises quickly to a plateau 
and then asymptotes to a constant as ^ oo. One can argue, along the lines of [22], that 
the delta-function spike at = had to be there because of the broken gauge invariance. 
But the existence of a gap is additional information, revealed by the calculations of 
but apparently not necessitated by symmetry principles. In BCS theory, the gap arises 
because of a pairing mechanism of otherwise nearly free quasi-particle excitations of a Fermi 
surface. No such mechanism is manifest in the gravity description; instead, the simplest 
way to characterize the gravity calculation is that photons with frequency less than 2A are 
very unlikely to penetrate through the condensate and be absorbed by the horizon. There is 
clearly something in common between BCS theory and the gravitational calculation, because 
the horizon represents the dynamics of the uncondensed charge carriers (i.e. the normal 
component), and absorption of a photon with > is associated with an excitation of these 
carriers. The obvious difference is that the charge carriers in the gravitational calculation 
(or, more precisely, the charge carriers in the appropriate holographic dual description) are 
strongly coupled even when they are in the normal state. 

The strong coupling inherent in gauge-string duals in the gravity approximation raises the 
appealing possibility that black hole constructions might provide useful physical analogies to 
the mysterious dynamics of electrons in high Tc materials that go beyond traditional ideas 
based on quasi-particle excitations of Fermi surfaces. But the black holes we study provide 
anything but a microscopic understanding of superconductivity: the gravity description is 
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more like Landau- Ginzburg theory, and the dual field theory would be the venue for some 
attempt at a microscopic theory comparable to BCS. 

Rather than presenting superconducting black holes as an incipient theory of high T^., 
we prefer the viewpoint that they are a new theoretical laboratory, seemingly divorced from 
traditional perturbative concepts, but capable of exhibiting assorted phenomena reminiscent 
of real materials. Perhaps a sufficiently comprehensive understanding of their dynamics will 
suggest new ideas which can also be applied successfully to real materials. 

The purpose of the present paper is to narrow the gap between black hole constructions 
and interesting high materials by introducing black holes with a p-wave gap. Although it 
is apparently a rf-wave gap that controls the dynamics of the cuprates, d-wave and p-wave 
are similar in that excitations of the normal component can be probed using low-frequency 
photons. 



3 The backgrounds 

We follow the conventions of [1] for the metric and gauge field, except for denoting the spatial 
directions of R^'^ as x and y rather than and x^. We will restrict attention to the limit 
of large g, where the metric is simply y4(iS'4-Schwarzschild, 



The gauge field ansatz is 



It is convenient to define 



1 - -| I rft^ + rfx^ + dy^ 



+ 



A = ^(r)T^dt + w(r)T^dx . 



dr'^ 



<i> = gL"^^ w = gL?w . 
If one also fixes a scale by setting r// = 1, then the relevant Yang- Mills equations are 

1 



(15) 

(16) 
(17) 



r 

^„ l + 2r3 ^, 
w H —w 



= 



r(r^ — 1] 

^2 



r{r^ — 1) (r^ — 1) 



where primes denote d/dr. These equations are similar to (B4) of [Ij because the ansatz 
f|T6|) is also similar. But in j^, the symmetry breaking term takes the form w{T^dx + r'^dy), 
which corresponds to wrapping the part of the gauge group generated by — call it f/(l)3 — 
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around the rotational symmetry group SO (2) that acts on x and y. Choosing instead 
w^T^dx — r'^dy) corresponds to wrapping (/{l)^ the other way around S0(2). We think of 
f|T6|l heuristically as a superposition of the two different wrapping solutions, in the same way 
that linearly polarized light is a superposition of left-handed and right-handed polarizations. 
This analogy has limited utility because the Yang-Mills equations governing the different 
"polarizations" are non-linear. 

In addition to breaking U{1)3, the condensate wr^dx picks out the x direction as special. 
Therefore, if back-reaction of the Yang-Mills field on the metric were included, then we would 
not expect to be able to set g^x = Qyy, as we did in (fT5l) . The wrapping condensate w{T^dx + 
r'^dy) is simpler in this regard, because although it breaks U{1)^ and 5*0(2) separately, 
it preserves a diagonal subgroup which makes the stress tensor isotropic in the x and y 
directions. 



The temperature of the horizon is 



47rL 



2 ' 



(19) 



where, as before, we have set th = ^- The total charge density p is proportional to the 
part of the electric field at infinity: if 

$ = Po + ^ + o(l), (20) 

then 

(") 

The charge density p„ in the normal component is proportional to the part of the electric 
field at the horizon: if 

$ = <l>i(r-l) + 0[(r-l)2], (22) 



then 



$1 

Pn = . (23) 

Far-field and near-horizon expansions for the rescaled fields $ and w take the form 

= ^0 + — + • • • 

Far field: { . ^ (24) 

w = h . . . 

r 
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$ = <|)^(r - 1) + . . . 
Near horizon: ^ (25) 
w = Wo + W2{r - 1)^ + . . . , 

and it is convenient to introduce rescaled versions of the total and normal component charge 
densities: 

~ 2 r2 

p = K gL p 



,.2„r2. Pi 



/ (26) 

We also define the superconducting charge density as = p — p„. A natural choice of 
order parameter is Wi, because the SU{2) currents dual to the gauge bosons have a 
symmetry-breaking expectation value 

rl s:a 



Jt) cx W^5l5l . (27) 



There is a one-parameter family of "allowed" solutions to the Yang-Mills equations (fT8|) . 
where allowed means that the far- field and near-horizon asymptotic forms, (124|) and (!25|) . 
are satisfied. Thus (fTSl) . ( |24l) . and ( |25l) specify a non-linear boundary value problem. To 
understand why there is only a one-parameter family of solutions, let us examine the far-field 
and near-horizon expansions separately. The generic solution to (ITS]) includes a constant 
term Wq in the far- field expansion of and this is disallowed because it corresponds to 
deforming the field theory lagrangian by some multiple of Jl- Another way to describe 
why Wq is disallowed is that if it is non-zero, then the condensate is not normalizable. In 
the expansion ([21]), all three parameters shown explicitly are independent, which matches 
a simple counting argument: four integration constants (for two second order differential 
equations) minus one (for the constraint Wq = 0) equals three. Requiring that the gauge 
field is smooth and well-defined at the horizon leads to the expansions (1251) . The parameters 
$1 and wq are independent, but W2 and higher coefficients can be determined in terms of 
them. Having only two independent parameters in the near-horizon expansion (i.e. $i and 
wq) means that there are two constraints at the horizon. Generically, these two constraints 
will be independent of the far-field constraint Wq = 0. So there are three constraints total on 
four integration constants, leading indeed to a one-parameter family of solutions. At special 
points, one of the horizon constraints may become degenerate with the far- field constraint, 
and this is when one finds two branches of solutions joining together. 

Solutions to the boundary value problem discussed in the previous paragraph can be 
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Figure 2: Each point along the contours plotted represents a solution to the non-linear 
boundary value problem specified by (fTSi) . (!24l) . and (125|) . Points on the line labeled "normal" 
are RNAdS solutions, and if charge density is held fixed, temperature rises as one moves to 
the left. Points on the curve labeled "superconducting" break the abelian gauge symmetry 
generated by f/(l)3. Points on the other curves also break the abelian gauge symmetry but 
are expected to be unstable. The point where the superconducting solutions join onto the 
normal solutions is labeled Tc because the simplest scenario is for there to be a second order 
phase transition at this point. 
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Figure 3: The fraction ps/p of the charge carried by the superconducting condensate and 
the order parameter Wi are plotted against the rescaled temperature Tj^fp. At Tc, Ps/p 
vanishes hnearly, while W\ vanishes as a/T^ — T. 



generated using a "shooting" procedure. First one guesses numerical values of $1 and wq. 
Then one uses the near-horizon expansion fl25|) to seed a finite-element differential equation 
solver, such as Mathematica's NDSolve. Next one matches the numerical solution to the 
far-field expansion ( l24l) . augmented by a constant term Wq. In this way one finds as 
a function of $1 and w^. The zeroes of this function correspond to the solutions of the 
boundary value problem: see figure [2l Hereafter we restrict attention to solutions where 
w(r) has no nodes. There are additional solutions with nodes, but one generally expects 
them to be thermodynamically disfavored, because radial oscillations in w cost energy. 

Thermodynamic quantities for solutions along the node-free symmetry-breaking branch 
labeled "superconducting" in figure [2] are plotted in figure [31 



4 Electromagnetic perturbations 

By making the black hole charged under the gauge symmetry generated by r^, we 

explicitly break the SU{2) gauge group down to U{1)^. We interpret U{1)^ as the gauge 
group of electromagnetism, which means that we plan to consider a weak gauging of this 
group in the boundary theory]^ As discussed in the introduction, the linear response to 
electromagnetic probes is described by the two-point function of the U{l)-i current, and in 
the dual black hole, this means that we want to know how linear perturbations of the 

^This situation is analogous to the Hubbard model, which has (at least) a global U{2) symmetry. The 
central U(\) is identified as the electromagnetic gauge symmetry, but electromagnetism is not explicitly part 
of the model. 
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component of the gauge field propagate. We persist in choosing the spatial momenta ki = 
in the x and y directions, so the photon is directed straight down into AdS^, as illustrated 
in figure [TJ 

As a warmup, we work out in section W7i\ the conductivity in two examples where it can 
be done analytically, including the normal state, where the condensate w is set to 0. In 
section 14.21 we explain how to set up the perturbation equations in the more difficult case 
of a symmetry-breaking background as described in section [3l In sections 14.31 and 14.41 we 
present results of a numerical study of axx{^) and o-yyiyj) which reveal a p-wave gap. 

4.1 Analytical calculations 

The simplest case to start with is pure AdS/^, corresponding to zero temperature, zero charge 
density, and no symmetry breaking. At the linearized level, the gauge coupling of SU{2) 
doesn't enter, so we will pass to free Maxwell theory in AdS4^: that is. 



S = ^ I d^x 



2k 



(28) 



The perturbation calculation is simple because the background geometry is conformally flat: 

ds"" = ^ {-df + dx" + dy^ + dz^) . (29) 

Equivalently, we may consider the line element (flSi) with th = 0: it is the same as (l29l) if 
one sets 

z = L^/r . (30) 

Conformal flatness is special because Maxwell's equations respect conformal symmetry. Thus 
the complexified photon perturbation is a plane wave: 

A^ = e-^'"(*-^) . (31) 



We chose the plane wave solution that moves in the positive z direction: that is, it moves 
away from the conformal boundary at 2; = and toward the degenerate Killing horizon of 
the Poincare patch, at 2 = 00. (In figure [H the positive z direction is downward.) This 
choice means that we will wind up computing the retarded Green's function rather than the 
advanced one. The Green's function can be read off from an expansion near the conformal 
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boundary: 

A, = e-*^*(l + iu;z + ...) = e-*"* (l + ^ + ...). (32) 



J, 

Comparing the last expression in fl32l) to f|T3|) . and using (jlj) and f|T^ . one finds 

<Txx = CToo = TT ■ (33) 

Because of rotation invariance, (jyy = a^x and a^y = 0. Hereafter we will normalize all 
conductivities against (Too by defining 

a., = ^ . (34) 

'-' oo 



Putting; (11), (til), (IMI), and (Bl) together, one has 



i Ax ^ 



L2 ^(0) ■ 



(35) 



A surprising result of [23] is that a^x = cfyy = 1 for the A(i5'4-Schwarzschild solution ( fTSjl . 
for all u and T. In the approximation where the gauge field (fT6|) doesn't back-react on the 
metric, this result persists so long as w = 0. The quickest way to derive it is to compute 
directly the linearized equation of motions for complexified gauge field perturbations of the 
background (fT5l)- (fT6l) : that is, A ^ A + a, where 

a = e-''^'al{r)T^dx . (36) 

The result of plugging this perturbation into the linearized Yang-Mills equations is 



2r^ + l ^ to^L^r^ 



< = , (37) 



r(r3-l) ' (r3-l)2j - 

where we have set r// = 1 as usual. Because the rotational symmetry is unbroken in the 
absence of the condensate, the same equation governs perturbations. The solution to flTTI) 
describing gauge bosons falling into the horizon at r = 1 is 

al = {r- l)"-/4-^(r2 + r + l)*-/«-^ | ^ ~^ ] , (38) 

2 2~ 



where we have used (1191) . The behavior oc (r — 1) ^'^/^'^^ is typical of solutions falling 
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into a finite-temperature horizon. The expansion of (138!) near the conformal boundary is the 
same as (1321) through order 1/r, so the conductivity is the same, as claimed. 



4.2 Electromagnetic perturbations of the superconducting phase 

In the presence of the condensate wr^dx, perturbations of the form (1361) mix with other 
components at the level of linearized equations. An ansatz which is sufficiently general to 
obtain consistent linearized equations is A ^ A + a, where 



(39) 



All the are functions of r. Plugging the perturbation (1391) into the linearized Yang-Mills 
equations, one finds that the mode obeys an equation of motion decoupled from the 
others: 



o2 Sr'^ 
d H 



-dr 



0. 



(40) 



1) ■ ■ (r3-l)2 r(r3-l) 

This is identical to (13) of [3], except that the last term has slightly different radial de- 
pendence. Unsurprisingly, the rescaled complex conductivity dyy exhibits similar gapped 
behavior to what was found in [3]: see figure HI Because the analysis is so similar to [3], we 
will not discuss it further here, but simply present the results in sections 14.31 and | 4.4[ 

The linearized Yang-Mills equations mix a^. with a\ and a^, resulting in three second 
order equations of motion. 



dl + -zdr + 



r(r3-l) ' (r3-l)2j - (r3 - 1 



2^2 



ai = 



^2"'^ 



d^r + -dr 



1 o 

a; + —at = 



r r{r^ — 1) 



r(r3 - 1) 

iujL?"w 
r{r^ — 1) 



(41) 



and two first-order constraints. 



■ r2/ 2\i 



iuL (a 



2/„l\/ 



2\i 



^'ai = 



II 



[wdr — w'] a 



0, 



(42) 



where, as before, primes mean d/dr. The constraints are not independent of the equations of 
motion: if one takes the r derivative of each constraint, the resulting second order equation 
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follows algebraically from the equations of motion and the undifferentiated constraints. It 
takes six constants of integration to specify a solution to the equations of motion, but two 
of them are used up in satisfying the constraints, leaving four independent solutions. Of 
these, two can be found in closed form and are related to residual gauge invariance, as we 
will discuss in more detail below. There is also a solution describing gauge bosons falling 
into the horizon, and another describing gauge bosons coming out. 

Let's focus on the infalling solution, which determines a retarded Green's function, as we 
have seen in easier examples above. Near the horizon. 



al = (r 



- 1)— [1 + ^3(1)^^ _ 1) _^ ^3(2) ^ ^ _ 



a] = (r - 1)— /4.^ ar\r - if + a^^r - 1) 
a2 = (r - 1)— 



(43) 



2(1)^ 1\ , 2(2)/ -,n2 

/ r — 1 + a/ r — 1 



and all the coefficients be determined once the background and u are specified. 

Near the conformal boundary, a generic solution to the equations of motion takes the form 



.3(1) 

3 ^srm , 



r 

Am 

al = + ^ + . . . (44) 
r 

42(1) 

a? = Af) + ^ + ..., 



and the constraints impose the relations 



• r2 I ~ 42(1) ~ a2(0) ^ 

zuL'Af^ - PoA]^'^ + p^Af^ + Af) = 



(45) 



where the coefficients ps and Wi are the ones appearing in fl24j) . The infalling solution is 
unique up to an overall scaling, which is fixed once we choose the leading behavior of 
to be (r - l)-»^/4^^ as in the first line of (H3|) . Thus the far-field coefficients Am^^ are in 
principle known as functions of u once the background is specified. We claim that 



n, ,7-2/12(0) , ~ 4l(0) 



The first term in parentheses is the expected result based on the considerations of 
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The second term has to do with solutions to (HT!) and (H2|) which are pure gauge outside the 
horizon, as we will now explain. 

An infinitesimal gauge transformation of the SU (2) gauge field takes the form 6 A = Da, 
where D = d + gA is the gauge- covariant derivative and a is an adjoint scalar gauge function. 
Let's consider the case 

a = e "'~"a"T 



After performing the split A ^ A + a oi the gauge field into background and fiuctuating 
parts, we can view the infinitesimal gauge transformation as acting only on a = e~*'^*a^r"(ix'^: 

5{e-'^'al) = a^(e-^"*a'^) + ^e'^^'^Aj^e-^^^a'^ . (48) 

If any a°' depends on r, then the gauge-transformed perturbations will include components a" 
which weren't present in the original ansatz fl5^ . Setting these components to zero amounts 
to choosing a form of axial gauge, and the gauge transformations that preserve axial gauge 
are the ones where a"" doesn't depend on r. Dependence on and is excluded because 
we are always considering modes with zero spatial momentum. We also set a'^ = because 
it would introduce components of the perturbations that are not present in the ansatz fl39p . 
To summarize: we are interested in infinitesimal gauge transformations of the form (1481) 
where and are constant and = 0. The explicit form of this restricted set of gauge 
transformations is 

Sal = '"^'5^ 

5al = -iujL^a^ - ^a^ (49) 
5al = -iuL'^dP' + <la^ , 

where in order to simplify notation we have defined a" = a"" /L'^. It is readily checked that 
the expressions in (149!) solve the equations of motion (14T!) and the constraints (1421) . This had 
to happen because (HT|) - (H2|) came from the gauge-invariant Yang-Mills equations. These 
are the two closed- form solutions which we mentioned in the text following (1421) . 

Up to an overall scaling, there is a unique linear combination of a^,, a\, and a^ which is 
invariant under the gauge transformation (H9l) : 

s ^iiuL'^af — ^aj 
0-1 = 4 + — ■ (50) 

The conductivity cFxx captures some gauge-invariant information about the bulk theory, and 
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as such it must be expressible in terms of a^. If one expands 

43(1) 

a^ = if) + ^ + ... (51) 
r 

near the conformal boundary, then the unique extension of (135!) that respects the gauge 
invariance is 

o'xx = — F?^wr- ^^'^^ 

UjU /l3(Uj 

This is precisely the result (1461) that we claimed earlier. In appendix [X] we describe how dxx 

fits into a 3 X 3 matrix of conductivities which can all be determined in terms of J^x^^ and 
43(1) 

4.3 Results of numerics 

Let us review the structure of the problem before discussing results. The gauge field back- 
ground f|T6l) is constructed by numerically solving the Yang-Mills equations f|T8l) in a fixed 
gravitational background, f|T5l) . subject to constraints near the conformal boundary and near 
the horizon, fl2^ and ^IM respectively. From a numerical solution, one can pick out the pa- 
rameters po, pi, and Wi appearing in (l46l) . A symmetry-breaking background with w > 
everywhere is labeled uniquely by the value of T/ ^/p, which has a maximum value of approx- 
imately 0.125. It is interesting that this value is within numerical error of 1/8, but we don't 
see any reason why it should be exactly 1/8. With a numerically constructed background in 
hand, one chooses a value of u;, initializes a finite-element differential equation solver close 
to the horizon using the series expansion (H3l) . solves fj4T]) . and extracts the coefficients A?x^\ 
^4^^^^ and A^-^^ appearing in ( l46l) by comparing the far-field behavior of the numeri- 

cal solution with the expansions (1441] . It is important to note that uj and L appear in the 
differential equations (14T!) and the conductivity formula (146!) only in the combination 

(53) 

where we have used ( IT9l) . (Recall that we have set rj^ = 1. If we had not, the left hand side 
of (!53l) would be instead cuL^/rj^, because then T = 3r///47rL^.) Thus it is more precise to 
say that one chooses a numerical value for the dimensionless quantity uo/T and determines 
CFxxi which is also dimensionless, in terms of it. One expects that for large enough cu/T, 
c^xx 1- This is because the condensate involves dynamics with a characteristic energy 
scale, which turns out to be ^/p. Provided we avoid the extreme limit T — > 0, i/p and T are 
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comparable. If ^ \fp, the propagation of the gauge bosons should be largely insensitive to 
the condensate: instead, its wave function approximately takes the form fl31l) that we found 
for photons in pure AdS^^, and a^^. ~ Ooo- 

Numerical computations can only detect the continuous part of axx{!^\ but there is also 
a distributional part with some interesting structure. Because a^xi^iS) is proportional to a 
retarded Green's function, it is analytic on the upper half-plane of complex lo. It therefore 
satisfies the Kramers-Kronig relations: 



The reason that a^x — 1 appears in (15^ rather than dxx itself is that it is Oxx — 1 which 
vanishes as oo, and such vanishing is a necessary condition in order to obtain fl54l) 

from a contour integral in the upper half-plane. V denotes the principle part of the integral. 
Evidently, a simple pole in lmaxx{^^) at u; = cuq implies a delta-function contribution 6{u!—uJo) 
to Keaxxi^^)- The positivity constraint on the real part of conductivities applies separately 
to the continuous and delta-function parts of Keaxxi^^), so any pole of Imaxxi^^) on the real 
axis must have positive residue. 

Plots of axx{^) and d-yyiu) are shown in figure HI The conspicuous features are: 

1. dxx and dyy both approach 1 as u; becomes large, as expected on general grounds. 

2. dyy displays gapped dependence similar to the findings of [3], with A ^ \\/p- That 
is. Re a is very small in the infrared, then rises quickly at = 2 A = uog -^/p, with 
a slight "bump" a little above Ug that is reminiscent of the behavior expected for 
fermionic pairing. We use the notation uog even though it's not clear that Hedyy is 
strictly zero ioi Q < uo < uOg. 

3. There is a pole in Ima^.^. at u; = cjq ~ 1.8-\/p. Its residue becomes small as one 
approaches T^. It's clear from (H6l) that this pole had to arise, with residue proportional 
to the order parameter Wi. it comes from the denominator of the second term, and 



As discussed following ( l54ll . there is a delta-function contribution to Yieaxx at a; = uiq, 
whose coefficient is proportional to the residue of this pole. This resonance is perfectly 




(54) 



(55) 
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Figure 4: Rescaled conductivities a^ja; and ayy as functions of frequency. The dotted curves 
are the best fits of the Drude model prediction fl56p to Kea^xi^) at low u. 
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stable even at finite temperature, but perhaps if we relax some of the limits we have 
taken (like large A^) it would acquire a width. 



4. Kea^x never goes as low as Keayy, and its rise toward 1 happens more gradually and 
at a somewhat larger value of c^;, on order Uq. 

5. The small u behavior of Re be parameterized very accurately in terms of the 
Drude model, which predicts 

(Jo 

ReClDrude = 5 (56) 

where Gq = ne^r/m is a constant related to the density of charge carriers, and r is the 
scattering time. 

We are especially interested in the low-frequency dependence of the conductivities. Our 
numerical results make it plausible but not certain that ayy is strictly zero below a finite 
value of uj when T = 0. However, neglecting the back-reaction of the gauge field may not 
be a valid approximation for very low temperatures. On the other hand, the narrow Drude 
peak in a^x suggests conductivity due to quasi-particles whose scattering time diverges as 
T — > 0. Putting the behavior of dxx and d-yy together suggests a very special type of "node 
in the gap," namely one which is infinitely narrow as a function of angle in Fourier spacej^ 

4.4 Fits of temperature-dependent quantities 

In order to extract some simple quantitative information from our numerical results, we 
considered the dependence of various dimensionless quantities on the rescaled temperature 



^Wc thank D. Husc and P. Ong for discussions that led to the picture of an infinitely narrow node 
described here. 
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Tj^fp. Our findings can be summarized briefly as follows: 



f« exp |o.303^ - 2.20 

Pn I T 




(57) 




lim ^= \m.axx{^) 0.52 




lim —=\m.a^AuS) ^ 0.55 



lim Im(Ta;j;(u;) 0.28 . 



a;— XiiQ 



The approximately equalities in (157|) are in some cases quite close over a substantial range 
of ^/p/T, and in others represent no more than a T — ^ extrapolation: see flgureO None 
of the relations flFTl) should be taken too seriously, because they were made over intervals 
where T / ^/p varied only by a factor of 5. A particularly challenging case is the quantity 
lim^^o Inio"^.^(a;). The u; ^ limit converges slowly because of a "shelf effect:" for values 
in a region around uj ~ 1/r, we observed -^Ima^x ~ 0.55 at low temperatures, which is 
the same value as we flnd in the u —>■ limit for Imayy. But for cu ^ l/50r, we observed 
instead the value some 6% smaller quoted in fl57p . Our numerical algorithms aren't optimized 
for extremely small T and u, and it's possible that this shelf effect goes away at very small 
T, so that the residues of Imaxx and Imdyy agree in this limit. But the balance of evidence 
from our numerical exploration is that this does not happen, or happens very slowly as T is 
decreased. 

5 Stability calculations 

We expected that the p-wave backgrounds (fT6!) would be unstable against small perturba- 
tions that would eventually turn them into backgrounds of the type studied in [1]. These 
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Figure 5: Temperature-dependent quantities and approximate fits, as explained in (1571) and 
the surrounding text. We have defined axx,o = hm^^o Reaxxli^), o'yy^ = lim^^oRe dyy^u), 
Res^=olm a^x/Vp = lim^^o ^Imaxxi^^), Res^=olmayy/^/p = lim^^Q ^lmayy{uj), and 
Res^=<^o Im axx/\/p = hm<^^^o Im axx{(^)- 



24 



backgrounds display behavior analogous to a p + ip gap But the opposite seems to be true: 
numerical explorations of quasinormal modes close to show that it is the p + ip-wave back- 
grounds that are unstable, and it seems that they evolve toward pure p-wave backgrounds, 
which are stable. In section [5?T] we exhibit the equations describing the perturbations of the 
pure p-wave backgrounds that we thought would be unstable and explain how the lowest- 
lying quasinormal modes exhibit stability instead, close to T^. In section 15.2^ we show that 
similar perturbations of the backgrounds studied in [1] exhibit an instability slightly below 
T 

c- 



5.1 Quasinormal frequencies of j9-wave backgrounds 

Let us begin by explaining why we thought p-wave backgrounds would be unstable. At 
T = Tc, both the r^dx mode and the r'^dy directions exhibit marginally stable modes. So a 
natural expectation is that both become unstable for T < Tc. Yet the p-wave backgrounds 
described in section |3] involve only r^dx, whereas the p + ip-wave backgrounds of [4] involve 
the combination r^dx + r'^dy. In the latter case we are taking advantage of both directions 
of instability, and it seems reasonable that such a configuration should be preferred. But 
this reasoning ignores the non-linearities of the Yang-Mills equations. It turns out that 
condensing in the T^dx direction stabilizes against condensation in the r'^dy direction — at 
least, close to T^. That stabilization is what we are going to address in this section. 

Starting from the backgrounds fITB]) . we want to study r'^dy perturbations, which is to 
say tty. At the linearized level, couples with a^, so we are forced to examine the combined 
perturbation A ^ A + a, where 

a = e-'"* {ay + ajr^) dy . (58) 

^"^The analogy to & p -\- ip gap is apt because the combination r^dx + r^dy distinguishes an orientation 
on and imphes a spontaneous magnetization. To see this, note first that the positive charge of the black 
hole under U{1)^ privileges over . The structure constants e"''^ of SU{2) then privilege the ordering 
(T^,r^) over (T^,r-'^), because having distinguished the positive direction in the Lie algebra lets us set 
c = 3. Finally, T^dx + r^dy "locks" this orientation in the Lie algebra to an orientation dxAdy on R^. More 
physically, a contribution w{T^dx + r'^dy) to A means that there is a term vP'T^dx A dy in F, representing a 
spontaneous magnetization that again picks out an orientation dxAdy in R^. In any case, the symmetries of 
this state are clearly those ot a p + ip gap whose ip component is of identical magnitude to its p component, 
so that the gap is uniform in magnitude but has a phase that rotates by 2tt as one goes once around the 
Fermi surface. 
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Re w/. 



Figure 6: Quasinormal frequencies corresponding to the perturbation ( 159!) of the p-wave 
superconducting background ( fT6|) near the critical temperature. The quasinormal mode 
spectrum is symmetric about the imaginary axis, and we are only showing the quasinormal 
frequencies with non-negative real parts. The arrows are in the direction of decreasing 
temperature, and the number displayed next to each quasinormal frequency represents T /T^. 
The blue points correspond to backgrounds with no condensate above Tc; the brown points 
correspond to backgrounds with no condensate below Tc; and the red points correspond 
to superconducting backgrounds below T^. The superconducting backgrounds also have a 
quasinormal mode at u; = (see main text) which is not displayed. The backgrounds with no 
condensate below have quasinormal frequencies with positive imaginary parts, indicating 
an instability. The other backgrounds (namely normal state above and superconducting 
below Tc) appear to be stable. 
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The equations of motion read 




2r^ + 1 



(r3- 1)2 



1 





a, 



2r3 + 1 



(r3 - 1)2 



(59) 





r(r3 — 1) 




0. 



The appropriate boundary conditions for quasinormal modes are that and should vanish 
at the boundary of AdS4^ and that a should be a function only of the infalling coordinate 
t + ^;^log(r — 1) at the black hole horizon (where, as usual, rn = !)• These conditions 
can be simultaneously satisfied only for certain complex quasinormal frequencies u. Since 
we assumed e"*"^* time dependence, quasinormal frequencies with negative imaginary parts 
correspond to stable modes, while those with positive imaginary parts correspond to unstable 
modes. Solutions with purely real u correspond to true normal modes of the system. From 
the symmetries of the equations (l59l) and of the boundary conditions described above, it 
follows that if a; is a quasinormal frequency, then so is —uj*. So let's restrict attention to 
quasinormal frequencies with Recj > 0. Figure [6] shows how the lowest-lying quasinormal 
frequencies behave as functions of temperature close to T^. Above T^, the normal state 
is stable, and the quasinormal modes come in degenerate pairs with the same u. As we 
mentioned earlier, there are two quasinormal modes that become marginally stable at T^: 
their frequencies go to zero. One of these modes, involving only a^, stays right at = 
below Tc on the superconducting branch. It is a Goldstone mode describing spatial rotations 
of the condensate. The other mode is stable on the superconducting branch below T^. What 
makes it stable is the ^^"3"^-^^ term in the second equation of fl59l) . This term is like a positive, 
r-dependent contribution to the mass term of the gauge boson. Dropping this term amounts 
to passing to the normal state below Tc, and our normal investigation showed that this state 

- 2 

is unstable. So the ^ "3"^,^^ term is the advertised stabilization mechanism, and it is evidently 
due to the non-linearities of the Yang-Mills equations of motion. 

5.2 Quasinormal frequencies of j9 + ip-wave backgrounds 

We now wish to show that the large gL limit of the p + ip backgrounds studied in [1] are 
unstable, at least for T close to Tc. The instability decreases the ip component of the gap 
and appears likely to lead the system into a p-wave state like (I16p . Our strategy is to find 
out what happens to the modes which are marginally stable at Tc as we go slightly away 
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from the critical temperature on the superconducting and normal branches. 
At large the gauge field ansatz for the circularly polarized backgrounds is 



A = $(r)r^dt + w{r) (r^dx + r'^dy) , 



(60) 



and it is again convenient to define 



w 



gL'^w. 



(61) 



In the large g limit there is no back- reaction on the metric, so the metric is simply ( llSp . The 
equations of motion for $ and w are similar to f|T8|) . They are given explicitly in (B4) of [1], 
and we will not reproduce them here. 

There are many ways in which one can perturb the background (160!) . but the perturbations 
that might show an instability towards converting p + ip into p should be of the form 



e-^'^*ai(rMx - r^dy) + e-''^'a2{T^dx + r^dy) . 



(62) 



The ai perturbation changes the relative magnitude of the p and ip components of the 
background ansatz (!60!) . Nothing in the ansatz (!60l) picks out whether r^dx or r'^dy is the 
p-wave part (as opposed to the ip part) so changing the relative size of these two components 
with a linear perturbation can be interpreted as decreasing the ip component without loss of 
generality. The a2 component is a 90° spatial rotation of the ai component. The linearized 
equations for ai and 02 are 



2r^ + 1 
r{r^ — 1] 

2r3 + 1 
r(r^ — 1) 



■dr 



$2 



+ 



+ 



1: 



ai 



0,2 



(r3 - 1)2 
2zu;L^r^l> 



0-2 = 



-ai = 



(63) 



The perturbations should take the form of infalling waves close to the horizon and should 
vanish at the boundary, as in the case of the linearly polarized backgrounds examined in 
the previous section. Only for discretely many quasinormal frequencies are these boundary 
conditions satisfied. 

When w = 0, equations ( l63l) are the same as the equations for and given in ( 159|) . 
so at zero condensate the quasinormal modes coincide with the ones displayed in figure O 
When T < Tc we find an instability whether or not there is a condensate: see figure [71 
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Figure 7: Quasinormal frequencies corresponding to the perturbation (162|) of the p + ip- 
wave background (160|) near the critical temperature. The spectrum of quasinormal modes is 
symmetric about the imaginary axis, and we are only showing the quasinormal frequencies 
with non-negative real parts. The arrows are in the direction of decreasing temperature, and 
the number displayed next to each quasinormal frequency represents T/Tc. The blue points 
correspond to backgrounds with no condensate above T^, the brown points correspond to 
backgrounds with no condensate below T^, and the red points correspond to superconduct- 
ing backgrounds below T^. The backgrounds with no condensate above T^, as well as the 
superconducting ones below T^, have quasinormal frequencies with positive imaginary parts, 
indicating an instability. The backgrounds with no condensate above Tc are likely to be 
stable. 
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This result could perhaps have been anticipated by noting that the w'^/r{r^ — 1) terms in 
(1^5]) enter with the opposite sign from the way they entered So instead of tending to 

stabilize perturbations, they tend to destabilize them. It's worth noting, however, that w is 
the coefficient of r^dx in fl59|) . whereas it is the coefficient of r^dx + r'^dy in fl63|) . So, the 
functional forms of w will differ in the two cases, becoming equal only in the limit T — T^.. 

6 Discussion 

The distinguishing feature of the superconducting black holes constructed in this paper is 
that the condensate is anisotropic, in the sense of picking out the x direction as preferred. 
This is in contrast to earlier constructions [U [21 [31 H] . What is special about the x directions 
is that the conductivity in this direction, a^x, becomes large at small but non-zero uo. So far, 
the situation is similar to p-wave superconductors. But in real materials, impurity scattering 
would keep axx finite for small non-zero cj, whereas in our setup, the only upper bound comes 
from the effects of finite temperature. The biggest difference from real materials — from the 
perspective of the electromagnetic response — is that ayy displays gapped dependence, similar 
to what was found for an s-wave construction in [3] . In real p-wave materials, the gap vanishes 
at ^ = and 6 = tt but has finite slope there. Gapped ayy suggests instead an infinitely 
narrow node in the gap: the slope of A as a function of 6 is infinite at 6* = and tt. To 
put it another way, the states which usually occupy a Dirac cone near a p-wave gap have 
been squeezed into a purely one-dimensional structure, at least in the limit of low energy. 
We emphasize that this picture of an infinitely narrow node in the gap is entirely heuristic, 
given that we do not have a microscopic description of the condensate in the language of a 
dual CFT. What we can say most clearly in the CFT language is that there is an SU{2) 
current algebra, and when there is a strong enough chemical potential for the charge density 
Jf, the component develops an expectation value. We are tempted to conjecture that 
•^m ~ i^ilm.Tij'ipi for some fermion fields ipi in a representation of SU (2). Then the condensate 
is composed of fermion pairs created by J^, which have one unit of angular momentum. 
Our results are preliminary in various ways: 

1. We didn't consider back- reaction of the gauge field on the metric. Back- reaction can 
be suppressed by taking the gauge coupling large, but this limit is non-uniform in that 
as T — > 0, the A], component of the gauge field gets larger and larger at the horizon, 
demanding a bigger value of the gauge coupling to justify the neglect of back-rection. 

2. Our conductivity calculations do not allow for spatial momentum. In other words, we 
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calculated a retarded two-point function Gr{ijJ, 0) of Jf at non-zero frequency but zero 
spatial momentum. A study of the electromagnetic response at non-zero k might help 
consolidate the heuristic Fermi-surface picture we have offered, or it might invalidate 
it and suggest a different interpretation. 

3. We encountered some curious numerical coincidences, ranging from Tc/^/p ~ 1/8 to 
the scaling of the "scattering rate" 1/r and the small uj limits of 1/ cixx and dyy approx- 
imately as rather than some fractional power of p„. The latter coincidence evokes 
the idea that the behavior of quantities like the scattering rate are largely controlled 
by kinematic factors of two incoming quasi-particles. It would be interesting if some 
of these numerical coincidences could be understood in terms of exact solutions to the 
Yang-Mills equations, or in terms of some systematic approximation scheme rather 
than brute-force numerics. 

4. The scope of our stability calculations is very restricted: not only have we limited 
ourselves to the no-back-reaction limit, but we also stayed close to T^.. Moreover, we 
do not claim to have considered every possible perturbation, only the ones that seemed 
obvious candidates for exhibiting instabilities. It would clearly be desirable to be more 
thorough. 

5. String theory or M-theory compactified on a manifold of positive curvature leads, often 
if not typically, to a theory of gauged supergravity which would include the Einstein- 
Yang-Mills lagrangian, with particular relations between the gauge coupling and the 
cosmological constants enforced by supersymmetry (assuming there is supersymme- 
try). Does superconductivity occur in any such construction? Is it the gauge bosons 
which condense first, or are there charged scalars which condense? Is there a higher- 
dimensional interpretation — for example, some sort of Gregory-Laflamme instability 
involving the compact dimensions? 

6. The constructions we have discussed probably generalize to higher dimensions. 

7. We have limited ourselves entirely to classical configurations, excluding any discussion 
of fluctuations. This would seem problematic in two spatial dimensions because of 
infrared divergences, but fluctuations are suppressed when the radius of AdS^^ is much 
larger than the Planck scale, corresponding to a large limit in the dual CFT. But to 
understand the condensate's contribution to the specific heat, presumably one should 
consider fluctuations. 
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We hope to report on some of these issues in the future. 
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A Conductivity and resistivity matrices 



The conductivities o^x and Oyy that we computed in section |4] are the cr^i^ and cr^^ entries of 
a generahzed conductivity matrix given by 

^;;nM = -G'^^'(^,0), (64) 

where 

G^'f(^,^) = -z 1 rf^xe-*-%(t)([j;^(t,f),4(0,0)]). (65) 

The matrix is block-diagonal, each block corresponding to a set of gauge perturbations 
in the bulk that mix with each other. To compute G^'n^ from the bulk, it is useful to pass 
to a complexified gauge perturbation, A — » A + a, and write the field strength perturbation 
as F — s> F + /. Then / = Da, where D = gA is the gauge-covariant derivative, and we 
need to work only to leading order in a. The quadratic action is 

-52 = ^ 1 rf^x (^-ICr''^ . (66) 

From now on let's focus on the 3x3 block corresponding to J^, J^, and — in other words, 
we keep only aj, of, and non-zero and set all the other to zero. Integrating fIBB]) by 
parts, we obtain 

S2 = ^j d^x [al* (eom for a^) + drj] , (67) 

where 

J = {aYd..a] + a-9.a?) + ^^j^<*d,al . (68) 

We henceforth set L = 1 to further simplify notation. 

As explained in section 14. 2[ we can find a solution to the equations of motion and two 
constraints for a\, af, and (that is, equations fHT]) and P2|) ) which is unique up to an 
overall multiplicative factor. The uniqueness arises because there's only one infalling solution 
at the horizon. But this purely infalling solution may not have the desired properties at the 
boundary of AdS, so we should consider gauge-equivalents of it: 

a\ a^a] — iua^ — $(5^ 

oPal - iua^ + <la^ (69) 
a"^ —* a a'^ + wa . 
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The "gauge parameter" a" is just a multiplicative rescaling factor, which we are free to 
include, and should be independent of r at least on some neighborhood of the boundary. 

In order to compute two-point functions of the Jj's, we want to prescribe the value of the 
gauge-transformed perturbations at the boundary by sending 



(70) 



a + wa 



01 = PI 



as r ^ CX3, and then take mixed derivatives d'^/d(3*d(3 of J evaluated at the boundary. 
Treating J\hAy as the only contribution to the on-shell action means neglecting possible 
contributions from the horizon, which can be justified along the lines of [T7j. To be precise: 
we start with a purely infalling solution a\, a^, and a^; then we send the boundary term 



through the gauge transformation (1691) with the a° chosen to satisfy (1701) : then we express 
the whole thing in terms of the /5^'s together with the far-field series expansion coefficients 
Am^^ of the original infalling solution. After all this processing, one finds 



lujcr 



(3, 



\Pil 



(71) 



with 



-11 -22 ^ PoPi 



a 



;;-23 



;;.3i 

^xt 



;;-32 



— CT- 



:i21 



Pi 



a. 



xt 



P^33 
(^xx 



i 



3(0) 



.1(1) ~ A Po + ^^^t 



pI-uj"^ 



.2(0)- . .1 

2(1) A 'PG-^^A^ 

.3(0) 1 A -I>l 
-t^x 



1(0) ■ 



i_^PoW]_ 
upl — u'^ 



pI-uj^ 



- AI^^'^po + iuA 



2(0)- 



-a;2 



(72) 
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a^l is the same as the one obtained in (H6|l because the calculation that led to (H6l) is a 
special case of the one above. The far-field limit of the two constraint equations fl45p enforce 



the relations al^ 



all and a?^ 



'tx 



-^ft, so 



Pi 



Pi 



PoWi \ 



Pi 



PoPi 



Wi 



PI 



- 

poWl 



upl 



Pi 



Wi 



\ 



upl 



Pi 



^33 

^ XX 



(73) 



with 5"^^ being given in (1721) . Since Pi^ and Wi are real, the only contribution to the 
hermitian (dissipative) part of cr comes from o"^^. 

It is worth noting that given the conductivity matrix cr one can compute its inverse 



After imposing we obtain 



Po 
;— 

Pi 



\ 



where p^^^ can be computed from 



533 

rxx 



Pi 


Pi 


-33^ 
rxx 




■ Po 





Pi 


Pi 




rxx 





Pxx J 



(74) 



i(i)i 



43(0) 



(75) 



Here l/pss is just the numerical reciprocal of P33, not a matrix inverse. Of the entries of p, 
pII. is the most interesting, because it represents E^./ with the constraint that = = 0. 
Using again the far- field constraints (H5|) . we obtain 



Im — 

Pxx 



Im a. 



-33 



W?Po 



ujpiipl - 



(76) 



Because Wi, po, and pi have no dependence, ( 176!) implies that Rel/p^j^ and Rea^^ differ 
only by a sum of three delta functions at u; = and ±po- The form of (!75l) suggests that 
Iml/pi^j^ doesn't have poles at u; = ±po- Numerical evaluations confirm this. So Rel/p^| 
doesn't have a delta function singularity at a; = ±po- 
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The singularity structure at a; = ±po is an unexpected feature of our calculations. Our 
heuristic understanding of it hinges on thinking of the boundary limit of as a chemical 
potential for the charge density J*^. The fluctuations that we are tracking in fITT]) include 
small rotations of the J*^ charge condensate into the J*^ and/or J*^ directions. Assuming 
that the total magnitude of the condensate remains fixed, this implies a slight decrease of the 
J*^ condensate. The chemical potential fights against this and tends to pull the condensate 
back to the J*^ direction, similar to gravity's pull on a pendulum. In the approximations we 
use, this resonance is undamped, but some kind of corrections — perhaps loop corrections on 
the gravity side — might damp it and smooth out the delta function singularity into a peak. 
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